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Here we report recent measurements on acetylene (C 2 H 2 ) ices at temperatures applicable to the outer 
Solar System and the interstellar medium. New near- and mid-infrared data, including optical constants 
(n, /<), absorption coefficients (a), and absolute band strengths (A), are presented for both amorphous and 
crystalline phases of C 2 H 2 that exist below 70 K. Comparisons are made to earlier work. Electronic ver- 
sions of the data are made available, as is a computer routine to use our reported n and k values to sim- 
ulate the observed IR spectra. Suggestions are given for the use of the data and a comparison to a 
spectrum of Makemake is made. 

Published by Elsevier Inc. 


1. Introduction 

The simpler hydrocarbons have long been known, or suspected, 
to exist in atmospheres and on surfaces in the outer Solar System 
where they have been studied remotely from Earth and Earth orbit 
and also by spacecraft instruments. Examples of hydrocarbons 
reported include frozen methane (CH 4 ) on Pluto, ethane (C 2 H 6 ) in 
cometary comae, Titan’s ethylene (C 2 H 4 ), and propane (C 3 H 8 ) in 
Saturn’s atmosphere. In addition to these four hydrocarbon mole- 
cules, acetylene (C 2 H 2 ) also has been observed in the outer Solar 
System. In cometary comae it is detected with an abundance rela- 
tive to H 2 0 of 0.1 -0.5%, and is thought to originate in the ices of 
cometary nuclei (Mumma et al., 2003). Such ices are believed to 
be representative of cold interstellar and outer protoplanetary disk 
regions where comets form. Trace amounts of gas-phase C 2 H 2 also 
are found in the hydrogen-rich atmospheres of all of the giant 
planets, and acetylene has been firmly identified as a major trace 
gas in the N 2 -dominated atmosphere of Titan. In fact, Titan’s atmo- 
spheric C 2 H 2 may participate in chemical pathways leading to 
polycyclic aromatic hydrocarbons (PAHs), compounds thought to 
play a role in the production of Titan’s aerosols. 

Acetylene’s astrochemistry also can include solid forms (ices) of 
the molecule on the cold surfaces of Pluto and other trans- 
neptunian objects (TNOs) where frozen C 2 H 2 may participate in 
chemistry driven by cosmic radiation. Laboratory work already 
has shown that C 2 H 2 and C 2 H 6 are radiation products of CH 4 
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(Bennett et al., 2006), and so both molecules could be present on 
outer Solar System icy surfaces with CH 4 . See Johnson (1989) for 
details of the radiation environment and time scales for such radio- 
lytic changes. Although the evidence for complex hydrocarbons on 
Pluto’s surface still is limited (Sasaki et al., 2005; Merlin et al., 
2010; DeMeo et al., 2010), frozen CH 4 has long been known to exist 
there (Cruikshank et al., 1976). Aside from Pluto, spectral features 
of CH 4 and C 2 H 6 ices are consistent with observations of Make- 
make (2005 FY 9 ) by Brown et al. (2007) and possibly of Quaoar 
too (Schaller and Brown, 2007). 

Beyond the Solar System, C 2 H 2 has been identified in the gas 
phase in molecular hot cores and protostellar environments, where 
its percent abundance relative to H 2 0 is 0.1 -0.5. This abundance 
cannot be matched by theoretical models without the inclusion 
of C 2 H 2 molecules subliming from icy grain mantles (Brown 
et al., 1988; Nguyen et al., 2002; Stauber et al., 2005). Carr and 
Najita (2011) have analyzed rovibrational emission spectra from 
the Spitzer Space Telescope to determine the gas-phase abundance 
of C 2 H 2 in circumstellar matter around T Tauri stars. See also 
Willacy and Woods (2009) for acetylene in protostellar disks. 

The acetylene molecule possesses a linear geometry with D^h 
symmetry and so it has no permanent dipole moment and there- 
fore no microwave spectrum. This means that remote sensing 
studies of extraterrestrial C 2 H 2 have relied heavily on infrared 
(IR) observations and lab-based work for their interpretation and 
quantification (de Bergh et al., 2008). Although IR studies of gas- 
phase C 2 H 2 have a long history (Levin and Meyer, 1928), laboratory 
work on C 2 H 2 ices is far more limited, which in turns hinders 
searches for solid-phase acetylene on TNOs. In early studies, 
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Bottger and Eggers (1964) measured the mid-IR spectra of acety- 
lene at 63 K and Dows (1966) reported a few measurements of 
the intensities of three IR features for acetylene ices at 65 K. The 
first paper to more completely quantify acetylene’s IR spectrum 
apparently was that of Khanna et al. (1988), which presented 
mid-IR measurements of crystalline C 2 H 2 at 70 K. Accompanying 
that work was a table of thirty-one n and k values constituting so- 
lid acetylene’s optical constants from 3250 to 710 cm -1 , and which 
remain all that are available in print. 

Other work on acetylene ices has appeared more recently, but 
none with direct measurements of band intensities. Boudin et al. 
(1998) published mid-IR spectra and peak positions for solid 
C 2 H 2 near 10 K, both neat and mixed with H 2 0 and CO ices. Infra- 
red band strengths were estimated by comparing acetylene band 
areas to the areas of IR features whose intensities were known 
from the work of others. In the near-IR region, publications on 
C 2 H 2 ices appear to be limited to a single spectrum covering 
7150-4000 cm -1 (Quirico et al., 1999). Most recently, Knez et al. 
(2012) reported mid-IR spectra of solid C 2 H 2 and C 2 H 2 mixed with 
H 2 0, CO, C0 2 , and CH 4 ices at multiple temperatures. From that pa- 
per, mid-IR spectra of pure amorphous C 2 H 2 made at 20 K and then 
warmed to 80 I< are available electronically. Knez et al. (2012) also 
measured positions, widths, and strengths for three mid-IR bands. 

Motivated by this scarcity of quantitative spectral data on solid 
C 2 H 2 and other ices, we recently have initiated a research program 
to measure optical constants of small molecules relevant to the 
chemistry of the outer Solar System. In the present paper we pres- 
ent a complete set of spectra and optical constants from roughly 
7000 to 500 cm -1 for both the amorphous and crystalline C 2 H 2 at 
multiple temperatures. In the next section we describe our exper- 
imental procedures for recording IR spectra and measuring the 
index of refraction of C 2 H 2 ices as a function of temperature. This 
is followed by details of our analysis to determine optical con- 
stants, band strengths, and absorption coefficients for amorphous 
and crystalline C 2 H 2 in the near- and mid-IR regions. Our labora- 
tory results then are presented, followed by optical-constants 
calculations and some possible applications. Finally, in Section 5 
we compare our near-IR spectra of C 2 H 2 with the spectrum of 
the dwarf planet Makemake (2005 FY 9 ). 

2. Laboratory procedures 

Our experimental setup, ice preparation method, and IR spectral 
techniques have been described (Moore et al., 2010), and so only a 
summary is needed here. Briefly, we used a Perkin Elmer Spectrum 
GX single-beam FTIR spectrometer set up for measurements from 
7000 to 500 cm -1 (1.43-20 pm). Over this entire range we used a 
KBr substrate and vacuum windows on a cryostat capable of cool- 
ing to 1 2 K. The optical constants of the substrate were taken from 
Li (1976). Spectra were measured at 2-cnrr 1 resolution using a 
mid-DTGS detector and a KBr beam splitter, and with 100 scans 
per spectrum. 

The acetylene used in our experiments was Matheson Research 
Grade, purified by distillation from a liquid nitrogen/ethanol slush 
bath (163 K) to remove acetone, used as a stabilizer, and any resid- 
ual H 2 0 and C0 2 . 

2.2. Ice preparation 

The same procedure was used in the near-IR and mid-IR regions 
for ice growth and spectral measurements. After cooling the blank 
KBr substrate, its spectrum was recorded at 12, 20, 30, 40, 50 60, 
and 70 K. Gas-phase C 2 H 2 then was condensed from a vacuum sys- 
tem onto the substrate to form an ice. Amorphous C 2 H 2 formed 
when the deposition was carried out at 12 K, and crystalline C 2 H 2 


when the condensation took place at 70 K. The resulting ice’s thick- 
ness (ft) was determined by laser interferometry using Eq. (1) 

h = Nm ngesA (1) 

where Nf ring es was the number of interference fringes observed with 
a laser of wavelength X passing through an ice of refractive index n x 
and with the laser beam’s path approximately perpendicular to the 
ice’s surface (Groner et al., 1973). In our work, X = 670 nm and so the 
corresponding refractive index is denoted n 670 . Ices grown for mid- 
IR spectra were limited to about two fringes (ft « 0.5 pm) to prevent 
the saturation of the stronger bands, whereas ices grown for near-IR 
spectra typically gave at least twenty fringes (~5 pm) so as to 
ensure good signal-to-noise ratios for the weaker near-IR features. 
Since we could estimate about 1 /8 of a complete fringe, the error 
in ice thickness was about 6% for the thinnest ices, dropping to 
about 0.6% for the thicker samples. The growth rate of C 2 H 2 ices 
was roughly 1 pm h -1 . 

2.2. Infrared spectral measurements 

Transmission spectra of amorphous ices were recorded from 
7000 to 500 cm -1 (1.43-20 pm) at 12 K, and then at 20, 30, 40, 
50, 60, and 70 I< as the sample was warmed. For crystalline ices, 
an initial spectrum was recorded at 70 K, and then at 60, 50, 40, 
30, 20, and 15 K as the sample was cooled. The warming or cooling 
rate was about 3Kmin _1 . Each ice’s single-beam spectrum was 
ratioed against the background (blank) spectrum of the KBr sub- 
strate recorded at the same temperature, maintained with a preci- 
sion of ±0.3 K. As discussed by Anderson et al. (1985) and Khanna 
et al. (1988), the IR and Raman spectra of C 2 H 2 ice formed at 70 K 
are suggestive of an orthorhombic crystalline structure. Although a 
transition to a cubic crystalline phase occurs at 133 K (Smith, 
1969), it was not studied in our work since C 2 H 2 ice rapidly sub- 
limes above 80 K in our vacuum system. 

2.3. Determination of n 670 for C 2 H 2 ices 

Values of n 670 for acetylene ices were needed both for ice-thick- 
ness measurements, with Eq. (1), and for the calculation of optical 
constants {vide infra). Our method for measuring n 670 of ices has 
been described (Moore et al., 2010), so again only a summary is 
given here. We employed a dedicated instrument separate from 
the set-up used for recording near- and mid-IR spectra for optical 
constants. Briefly, we condensed C 2 H 2 at a constant rate onto a 
cooled KBr substrate within a vacuum chamber while simulta- 
neously, during the deposition, recording the interference fringes 
generated from two 670-nm lasers’ beams reflected at angles of 
0i = 6.758° and 0 2 = 44.04° to the substrate. The ratio, G/t 2 , of the 
oscillation periods of the two lasers’ reflected intensities gave the 
refractive index through Eq. (2). 

sin 2 (h - (fi/f 2 ) 2 sin 2 0, 

l-(t,/t 2 ) 2 

See, for example, Tempelmeyer and Mills (1968) for a similar use of 
this same equation. For each ice sample, the maxima and minima of 
more than ten fringes were measured for each laser so as to reduce 
the error in determining the periods of oscillation. The constancy of 
the acetylene’s deposition rate was supported by the high correla- 
tion coefficients (>0.99) of plots of the time (appearance) for each 
fringe extrema against each extrema’s ordinal number. See also 
Romanescu et al. (2010) for similar work. A Mattson FTIR spectrom- 
eter was interfaced to our two-laser system to verify the amorphous 
or crystalline nature of each ice. 
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The next section shows that to calculate optical constants a ref- 
erence value, n 0 , for the index of refraction of each acetylene ice 
was required. Since a slight reduction in n 670 was expected when 
extrapolating to longer wavelengths, we took n 0 to be 0.99 n 6 7 o 
based on a similar correction by Edwards and Philipp (1985) and 
our own work (Moore et al., 2010). 


3. Spectral intensities and optical constants 

For some icy solids, the positions, widths, and shapes of spectral 
features are sufficient for identifying molecules, particularly when 
coupled with some knowledge of the chemistry of the ice, such as 
the behavior of similar systems or isoelectronic species. However, 
the quantitative spectroscopic analysis of ices requires additional 
information and is usually more difficult to accomplish. At least 
three approaches to quantifying IR spectral features are in the lit- 
erature, namely absorption coefficients, band strengths, and opti- 
cal constants, but seldom are these found in a single study, such 
as the present one. 

Before outlining our work on IR intensities, we point out that 
spectroscopists have long recognized that optical complexities, 
both instrumental ones (Wilson and Wells, 1946; Ramsay, 1952) 
and those inherent to thin samples (Maeda and Schatz, 1961), 
can complicate band-strength measurements. In what follows, 
the adjectives “absolute” and “apparent” designate intensities that 
have and have not been, respectively, corrected for such matters. 
The use of these adjectives will perhaps seem reasonable to those 
in astronomical fields, but while these same expressions have been 
employed by some IR workers (e.g., Wexler, 1967; Scrocco et al., 
1972), others have preferred a slightly different terminology (e.g., 
Thomas and Orville-Thomas, 1971). Similarly, the notation we 
use for absolute and apparent IR intensities is the same as that of 
Maeda and Schatz (1961), and resembles that already in the astro- 
chemical literature, but some workers have preferred different 
symbols (e.g., Avram and Mateescu, 1972; Overend, 1982; Paul 
and Ford, 1989), and there appears to be no widely-accepted IUPAC 
convention to cover all such cases. 


3.1. Absorption coefficients and band strengths 

To a first approximation, the intensity of a beam of light before 
(J 0 ) and after (J) passing through a sample of thickness h can be 
described by I = I 0 exp(-a'fz), the Beer-Lambert law, where oc' is 
called the apparent absorption coefficient of the sample. 
Rearrangement and modification of this equation gives logi 0 (Jo/ 
J) = a'h/ln( 1 0), and so it follows that from measurements of both 
logio(/o/0 and Jz, a' can quickly be calculated, such as by taking 
the slope of a graph of logi 0 (J 0 /J) as a function of thickness h. Since 
logio(Jo/J) varies with wavenumber, due to the presence of spectral 
features, an absorption coefficient a' is a measure of spectral inten- 
sity at a particular wavenumber v. However, since IR bands extend 
over a range of wavenumbers, it is sometimes both necessary and 
more valuable to discuss spectral intensity in terms of the area 
under an absorption feature, specifically / a'(v)dv, which is the 
same as ln(10) /log 10 (J 0 /J)dv/Jz. Taking p N as the concentration 
(molecules cm -3 ) of absorbers in a sample then / a'(v)dv = p N A\ 
where A' is a measure of the band’s strength on a per molecule ba- 
sis and usually just termed an apparent “A value”, here in units of 
cm molecule -1 . All of these terms can be combined to give Eq. (3) 



for determining intensities of IR bands from spectra. The integral 
itself is simply the area of a spectral band recorded on a 


common-logarithmic scale. See Hollenberg and Dows (1961) for 
additional details and examples. 


3.2. Optical constants 


Although a' and A' are readily found by measuring peak posi- 
tions and band areas, these measures of spectral intensity include 
losses due to reflection and interference effects, which can be sub- 
stantial for thin ice samples (Maeda and Schatz, 1961). In this pa- 
per we designate absorption coefficients and band strengths that 
are free of such complications as a and A, respectively. These quan- 
tities can be found through a determination of the real and imagi- 
nary parts of a sample’s complex index of refraction m = n - ik 
where n and k together are known as optical constants. The quan- 
tity k is related to an absorption coefficient by a = 47iv/<. Values of 
both n and k depend on v and can be used to generate synthetic 
spectra with which to compare to astronomical observations. 

For the present paper we used an iterative process to determine 
optical constants from high-quality transmission spectra of C 2 H 2 
ices recorded in the near-IR (~7000-3500 cm -1 ) and mid-IR 
(4000-500 cm -1 ) regions for amorphous and crystalline samples 
at multiple temperatures. As already mentioned, the substantial 
differences in intensities between the fundamental (mid-IR) and 
higher- wavenumber (near-IR) features required that ices with 
two different thicknesses, “thin” for mid-IR and “thick” for near- 
IR, be grown for all amorphous and all crystalline C 2 H 2 spectral 
studies. Each single-beam spectrum, J(v), of an acetylene ice was 
divided by a single-beam spectrum, J 0 (v), of the blank KBr sub- 
strate to remove instrumental characteristics before conversion 
to an absorbance spectrum (i.e., -log 10 [J(v)/J 0 (v)]). As explained 
in Moore et al. (2010), the baseline of each such spectrum then 
was straightened, flattened, and brought to zero, if necessary, then 
converted back to an J(v)/J 0 (v) ratio, and then multiplied by a the- 
oretical channel-fringe pattern calculated from the thickness and 
n 670 of the ice. The result of these manipulations was a starting 
spectrum T(v) for a Kramers-Kronig analysis to obtain infrared 
n(v) and k(v) values. This is essentially the procedure that Hudgins 
et al. (1993, Section 3.2) used to correct for scattering and non-lin- 
ear baselines. 

Our iterative calculations of optical constants began with the 
lab spectrum T(v) just described, from which estimates of a(v) 
for the ice sample were made with 

a(v) = -llnT(v) (4) 


followed by the calculation of fc(v) from 


m 


«(v) 

4nv 


( 5 ) 


and then the first estimate for n(v) from Eq. (6), the Kramers-Kronig 
transformation, using the previously-determined n 0 = (0.99)n 6 7o. 


n(v) = n 0 + 



v'fe(v') 

(vf - (V ) 2 


dv' 


( 6 ) 


These n(v) values and the original k(v) values of (5) were combined 
to calculate transmission and reflection coefficients (e.g., see 
Hudgins et al., 1993) to then compute a synthetic spectrum, 
J(v)/J 0 (v). The latter gave a new set of a(v) and, from Eq. (5), a 
revised set of k(v) values, which through Eq. (6) gave a revised set 
of n(v). After each such cycle 

spectrum — > a(v) — ► k(v) — ► n(v) — > coefficients -► spectrum 


the new transmission spectrum was compared to the original spec- 
trum, T(v). The ratio of the calculated and original reference spectra 
was brought closer to unity by iterating until a self-consistent set of 
n(v) and fc(v) values was reached and the calculated-to-measured 
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ratio of spectra deviated from 1 by less than 10 -5 at all v values. To 
reach this level of agreement, as many as one hundred iterations 
were sometimes required for mid-IR spectra, but only twenty for 
the near-IR data. See Moore et al. (2010) for more on the spectral 
manipulations and the corresponding calculations, and Keefe et al. 
(2012) for details and references on the development and applica- 
tion of such methods. All computations were performed with Math- 
Cad software and selective independent checking with IDL routines. 

Having determined fc(v), we used Eq. (5) to calculate and plot 
a(v) as a function of v to give a final computed IR spectrum. Inte- 
gration of the IR bands of this spectrum, and application of Eq. 
(7), a modification of Eq. (3), then gave us absolute band strengths 
directly from our optical constants. 

A = f [ a(v)dv (7) 

HN Jband 

Again we note that these absolute band strengths (A) are different 
from the apparent ones (A') of Section 3.1. 

The largest source of uncertainty in A is the ice’s density. We as- 
sumed a mass density of p = 0.76 g cm -3 from work on orthorho- 
mic C 2 H 2 at 131 K by McMullan et al. (1992). Data on C 2 H 2 at 
lower temperatures is unavailable, but measurements on C 2 D 2 
down to 1 5 K by the same research group suggest that the density 
of C 2 H 2 will change by less than 10%. The concentration (number 
density, p N ) needed for Eq. (3) is readily found from p by p N = pN A l 
(MW), where N A is Avogadro’s number and MW is the gram molar 
mass of acetylene (26.04 g mol -1 ). Eq. (7) shows that if new p N val- 
ues become available in the future then they can readily be used to 
rescale our A values. 

4. Results 

We first present IR spectra of amorphous and crystalline acety- 
lene, then we examine the influence of temperature on the near- 
and mid-IR absorbances of acetylene, and then we present the re- 
sults of our n 670 measurements for C 2 H 2 ices. 

4.1. Infrared spectra for C 2 H 2 ices at 7000-500 cm 1 

Fig. la and b shows the IR spectra of C 2 H 2 ices at 12 and 15 K in 
the near- and mid-IR regions, with peaks and assignments listed in 
Table 1. These spectra quickly illustrate the differences in intensity 
of features in the near- and mid-IR regions. For both solid phases, 
the strongest band is the v 3 vibrational mode near 3230 cm -1 , 
which approached saturation for ice thicknesses near 0.5 pm. Band 
strengths and absorption coefficients for the v 3 feature and others 


are presented later in Section 5.2. We note here that our identifica- 
tions of the near-IR overtone and combination bands are taken 
from the gas-phase absorptions close to the positions we observed 
in the solid phase. More work is required for firmer assignments of 
these weaker ice features. 

Fig. 2 shows spectra of three near-IR and two mid-IR bands on 
expanded scales. In each case we compare the amorphous ice made 
at 12 K with a crystalline ice grown at 70 K and cooled to 15 K. The 
near-IR spectra of Fig. 2a, b, e, and f are for ices with a thickness 
near 5 pm, while the ices giving the mid-IR spectra of Fig. 2c, d, 
g, and h had a thickness near 0.5 pm. Two spectral properties of 
amorphous C 2 H 2 ice are readily illustrated in the mid-IR spectra. 
First is the symmetric shape of the single-peaked v 3 band, which 
has a maximum about 14 cm -1 greater than that of the peak for 
the crystalline ice. Second, the position of the peak of the strongly 
asymmetric v 5 band is about 28 cm -1 less than that of the crystal- 
line ice’s double-peaked v 5 feature. This pronounced asymmetry of 
an intense band at the smaller wavenumber (longer wavelength) in 
amorphous C 2 H 2 is reminiscent of a similar observation by Bertie 
et al. (1994, 1995). 

The spectra of Figs. 1 and 2 were checked with acetylene ices of 
various thicknesses to ensure that the features shown did not suf- 
fer from saturation effects in the regions of interest. That they did 
not was suggested by the linearity of the associated Beer’s Law 
plots (band area vs. ice thickness). We also recorded spectra at 
both 2.0 and 0.5 cm -1 resolutions to be sure that these same IR 
bands were not resolution limited. Only in the case of the sharp 
v 3 feature was a difference in area seen, and even then it was only 
on the order of 3%. 

In the early stages of our work two complications arose that 
should be noted. First, we found that traces of H 2 0 vapor in our 
vacuum systems were manifested by a pair of small IR features 
of the H 2 0-C 2 H 2 complex at 3642 and 3565 cm -1 (2.746 and 
2.805 pm). In the worse case, these H 2 0 molecules amounted to 
as much as 1% of the sample, but fortunately this could be reduced 
to the level of our spectral noise by prolonged pumping on and 
baking of our vacuum chambers. At the level of such noise, the 
upper limit to any H 2 0 present would have been on the order of 
0.001% of the acetylene sample. Engdahl and Nelander (1983) also 
detected a H 2 0-C 2 H 2 complex even in condensed mixtures of C 2 H 2 
diluted and trapped in dry research-grade argon. A second early 
observation was that when making C 2 H 2 ices at the lowest temper- 
atures too high a deposition rate gave spectral features that resem- 
bled those of our annealed ices. Reproducible IR spectra for 
amorphous C 2 H 2 , such as those shown as in Figs. 1 and 2, required 
deposition rates no higher than what were stated in Section 2.1. 



Wavenumber (cm 1 ) 


1 


V 4 + V 5 


Arrorphcus C 2 H 2 at 12 K 


Jj 


0.5 


Crystalline C 2 H 2 at 15 K 


A. 


A 


4000 3500 3000 2500 2000 1500 1000 500 

Wavenumber (cm 1 ) 


Fig. l. Near and mid-IR spectra of amorphous and crystalline C 2 H 2 at 12 and 15 K, respectively. 
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Table 1 

Peak positions of IR absorptions of solid C 2 H 2 . 


Amorphous, 1 2 K 

v (cm -1 ) 

Crystalline 1 5 K 

v (cm -1 ) 

Vibrational mode 3 

6465 

6485, 6482, 6472, 6466, 6454, 6452, 6442, 6431, 6421 

Vi + v 3 (v 2 + v 3 + 2v 4 ?) 

- 

5275 w 

Vi + 2v 4 + v 5 

5189 w 

5183 sh, 5174 

5v 4 + 3 v 5 

- 

4684 w 

Vi + v 4 + v 5 

- 

4502 w 

v 3 + 2v 4 

4072 

4077 

Vi + v 5 

3948 

- 

2v 2 

3862 

3864 

v 2 + 2v 4 + v 5 

- 

3438, 3428, 3393 sh 

v 2 + 2v 5 

3328 

3349 sh, 3333, 3318 sh 

Vi 

- 

3306, 3296 sh 

- 

- 

3279 

- 

3240 

3226 

V 3 

1419 sh, 1371 

1419 sh, 1393, 1374 sh 

v 4 + v 5 

742 

792 sh, 770, 761, 745 sh 

V 5 

634 

- 

V 4 


a Assignments are based largely on gas-phase data in the GEISA 2011 spectroscopic database (http://ether.ipsl.jussieu.fr/etherTypo/?id=1293&L=0, Laboratoire de Mete- 
orologie Dynamique, Ecole Polytechnique, France, accessed 29 May 2013) and in Herman et al. (2003), with assistance from Plyler et al. (1963), Bottger and Eggers (1964), and 
Anderson et al. (1985). Some assignments, particularly those for v > 4000 cm -1 , should be regarded as tentative. 







Fig. 2. Details of several major IR features of amorphous (a-d) and crystalline (e-h) C 2 H 2 at 12 and 15 K, respectively. 


4.2. Infrared spectra of C 2 H 2 ices as a function of temperature 

The near- and mid-IR spectra of C 2 H 2 ice on warming from 12 K 
are shown in Fig. 3. Between 40 and 50 K irreversible changes in 
band shapes and positions, particularly in the 6500, 3300, and 
750 cm -1 regions, suggest that the amorphous ice converted into 
a more crystalline form. Fig. 4 is an enlargement of the 
3300 cm -1 (v 3 ) feature clearly showing the change in the spectra 
of amorphous C 2 H 2 on warming from 40 to 50 K. However, the 
spectra of such warmed amorphous C 2 H 2 ices never matched the 
spectra of crystalline C 2 H 2 ices grown at 70 K, and this is illustrated 
in more detail in Figs. 5 and 6. The three panels of Fig. 5 show the 
v 5 band of C 2 H 2 ice made under different conditions. Fig. 5a shows 
IR spectra of amorphous C 2 H 2 made at 12 K and warmed to 70 K. At 
12 I< the peak of the v 5 band is at 742 cm -1 , but on warming it 
broadens and develops several components with a maximum at 
758 cm -1 at 70 K. Fig. 5b shows the spectrum of crystalline C 2 H 2 


formed at 70 K and then cooled. The spectrum of the 70 K ice has 
a double-peaked shape with maxima at 769 and 761 cm -1 . Only 
minor changes in intensity, shape, and position are seen on cooling 
the sample. Fig. 5c compares spectra of ices grown and measured 
at temperatures from 12 to 70 K, showing that each possesses a un- 
ique shape in the 20-40 K region. 

The contrast between the shapes and positions of the v 5 feature 
at 20 K using three different formation conditions is shown in 
Fig. 6. Spectrum (a) is for an ice grown at 12 K and warmed to 
20 K, spectrum (b) is for one grown at 12 K and warmed at 70 K be- 
fore recooling to 20 K, and spectrum (c) is for an acetylene ice that 
was grown at 70 K and then cooled to 20 K. The differences in (a), 
(b), and (c) are obvious. For acetylene we found from the v 5 feature 
in Fig. 6, and from every major band at 7000-500 cm -1 , that we 
were unable to form a purely-crystalline phase simply by anneal- 
ing an amorphous ice. We suspect that in each attempt a fraction 
of amorphous material remained. The purely-crystalline phase 
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Fig. 4. The v 3 C 2 H 2 absorption for amorphous C 2 H 2 formed at 12 K and warmed to 
70 K. 

required a higher formation temperature (i.e., 70 K) of the ice, 
which is why we chose to make depositions at 70 K to generate 
crystalline acetylene. 

4.3. Refractive indices for C 2 H 2 ices at 670 nm 

Fig. 7 shows the values we measured for n 670 for C 2 H 2 ices 
grown at 11-70 K, with each point being the average of at least 
three measurements. Using a standard method for the propagation 
of errors (Taylor, 1997; Garland et al., 2003), we estimated the 
uncertainties of these values to be about ±0.015, which is ~1% of 

^ 670 - 

The most-striking thing about Fig. 7 is undoubtedly the dip at 
20-40 K. Infrared spectra of the ices used for our n 670 measure- 
ments were recorded, and they had the same features already de- 
scribed, namely those in Fig. 5c in which C 2 H 2 was deposited at a 
similar set of temperatures. In other words, the fall and rise in 
n 670 on going from 20 to 40 K in Fig. 7 was correlated with the tem- 
perature-dependent spectral changes in Fig. 5c, where unique an- 
nealed forms of the acetylene ice were formed on making the 
amorphous material. At present we have no firm explanation for 
this behavior in n 670 , although it is possible to invoke the formation 
of inhomogeneous regions or even voids in the ice on formation. 
Other ices examined (e.g., N 2 , C0 2 ) did not show this same 
behavior. 





Fig. 5. The band shape and position of the v 5 C 2 H 2 absorption for different 
formation conditions: (a) amorphous C 2 H 2 formed at 12 K and warmed to 70 K, (b) 
crystalline C 2 H 2 formed at 70 K and cooled to 12 K, and (c) C 2 H 2 ice formed at 
different temperatures from 12 to 70 K. The thickness of each ice was about 0.5 pm. 
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Fig. 6. Different formation conditions result in different v 5 band shapes for C 2 H 2 
ice: (a) amorphous C 2 H 2 at 15 K, (b) amorphous C 2 H 2 warmed to 70 K and then 
returned to 15 K, and (c) crystalline C 2 H 2 cooled to 15 K. 



Fig. 7. Index of refraction of solid C 2 H 2 at 670 nm as a function of the ice’s growth 
temperature. 

In our calculations of ice thicknesses and in our optical con- 
stants work we were concerned only with ices grown at tempera- 
tures on the opposite ends of Fig. 7, highly amorphous C 2 H 2 formed 
at 12 K and highly crystalline C 2 H 2 made at 70 K. For this reason we 
used only n 670 = 1.33 and 1.36 in all calculations, the choice 
depending on whether the sample was deposited at 12 or 70 K, 
respectively. Our own estimates showed that varying the initial 
value of n 670 in a Kramers-Kronig analysis by 10% could lead to 
absolute errors in the final values of a and A of up to 1 0%, but the 
precise direction and amount of variation were difficult to predict. 

5. Discussion 

5A.IR spectra of C 2 H 2 ices at 7000-500 cm _1 

Our spectra of crystalline C 2 H 2 ices in Figs 1-6 are essentially 
the same as those published by others (Bottger and Eggers, 1964; 
Khanna et al., 1988), although our work has the advantage of cov- 
ering wider ranges of temperature and wavenumber. For amor- 
phous C 2 H 2 , the only published mid-IR spectrum of which we are 
aware is that of Boudin et al. (1998). Fig. 8 compares spectra of 
our (a) amorphous and (b) crystalline ices with a spectrum of 
Boudin et al. (1998) said to be for a sample of acetylene deposited 


at 10 K. The close resemblance of traces (b) and (c) in Fig. 8 sug- 
gests that it is likely that the ices of Boudin et al. (1998) were an- 
nealed, perhaps unintentionally, or that their C 2 H 2 deposition rate 
was high enough to have warmed or even partially crystallized 
their sample. This conclusion is supported by an inspection of 
the relative intensities and positions of their v 3 and v 5 absorbances 
compared to what we show in this paper. 


5.2. Optical constants and band strengths of C 2 H 2 ices at 7000-500 
cmr 1 

The IR spectra and other data presented in Section 4 allowed us 
to calculate optical constants of our acetylene ices. Figs. 9 and 10 
show n(v) and k(v) at seven temperatures for amorphous and crys- 
talline C 2 H 2 . Near- and mid-IR transmission spectra also are shown 
for each phase, and are included to illustrate the relative intensities 
of the molecule’s absorptions. These are our spectra that have been 
baseline corrected, using the steps outlined in Section 3, to give a 
maximum transmission of 1.00 over the range shown. For the 
amorphous ice, the sample for near-IR measurements had a thick- 
ness of 6.5 pm, while the sample for mid-IR spectra had a thickness 
of 0.50 pm. A value of n 670 = 1.33 was used for determining these 
thicknesses and in the subsequent Kramers-Kronig analysis at 
each temperature. For the crystalline ice, the sample for near-IR 
measurements had a thickness of 4.6 pm, while the sample for 
mid-IR spectra had a thickness of 0.49 pm. A value of n 670 = 1.36 
was used for determining these thicknesses and in the subsequent 
Kramers-Kronig analysis at each temperature of the crystalline 
material. Calculations were done at one point per 0.5 cm -1 . 

Table 2 summarizes some key results from our spectra and cal- 
culations, listing peak positions, absorption coefficients (a), and 
band strengths (A) for C 2 H 2 ices near 12 K. An analysis of such val- 
ues from 12 to 70 K showed that neither a nor A possessed a strong 
temperature dependence (<10% change) for either the amorphous 
or crystalline phases of C 2 H 2 ice. For example, cooling from 70 to 
15 K gave little change in A for the v 5 band of Fig. 5b, and neither 
did warming an amorphous sample from 12 to 40 K, as in Fig. 5a. 
However, on warming an amorphous sample through the 40-to- 
50 K region A(v 5 ) increased by about 22% when the sample crystal- 
lized. In short, the a and A numbers of Table 2 can be used for a first 
approximation of a and A at other temperatures provided that 
there is no phase difference involved. 



Fig. 8. Spectra of (a) amorphous and (b) crystalline C 2 H 2 at ~12 K from this work 
compared to (c) a spectrum of C 2 H 2 reported to be deposited and measured near 
10 K (Boudin et al., 1998). 
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Fig. 9. (a) Near-IR (left) and (b) mid-IR (right) transmission spectra and optical constants of amorphous C 2 H 2 . The vertical bar indicates the scale, and spectra are stacked for 
clarity. The samples giving the near- and mid-IR features shown had thicknesses of 6.5 and 0.50 pm, respectively. 


Another comparison we made was between A and A', the abso- 
lute band strengths determining by integrating k values and the 
apparent band strengths calculated by integrating the laboratory 
spectra. Table 3 compares our values of A and A for amorphous 
acetylene, and shows that the differences are small. As already sta- 
ted, we know of no comparable results for amorphous C 2 H 2 in the 
literature, but for crystalline C 2 H 2 Table 4 lists our A and A values 
alongside those from two other groups. Our results are in reason- 
able agreement with the work of Khanna et al. (1988), but little 
more can be said since the temperatures they used to make their 
ices are unknown, as is the nature of the annealing treatment they 
used. Dows (1966) did not explicitly state either the spectral reso- 
lution of his measurement, his sample’s thickness, or the substrate 
material. His A'(v 3 ) is nearly twice the size of ours, and Khanna’s, 
but other values are comparable. Beyond these general statements, 
little can be said in a quantitative way. Checks for possible satura- 
tion effects and resolution limitations were not mentioned in 
either of the previous papers, although we presume that such 
checks were made. 


An important test of our optical constants is their ability to pre- 
dict spectral variations among ices of different thicknesses, as op- 
posed to just those particular thicknesses used for our Figs. 1-6. In 
most cases such differences were small, and amounted to the usual 
increases in intensity with thickness. These were easily replicated 
with the usual J = J 0 exp(-4nvkh) Beer’s Law relation using our 
computed k(v) values. However, the v 5 feature of crystalline acet- 
ylene proved to be a more-severe challenge. Fig. 11a compares 
spectra of crystalline C 2 H 2 ices calculated with our fc(v) values 
and Beer’s Law at six thicknesses. The changes shown are as ex- 
pected, but these spectra do not match the laboratory data. In con- 
trast, Fig. lib shows IR spectra calculated with equation Al of 
Swanepoel (1983), which uses both n(v) and k(v) to account for 
the optical interference and reflection effects present. Significant 
changes in the band’s shape and the peaks’ positions are obvious 
on going from 0.25 to 2.00 pm. We emphasize that it is the calcu- 
lated spectra of Fig. lib that match our laboratory spectra to with- 
in 10 -5 , as mentioned earlier. Only when both n(v) and k(v) were 
used to correct for optical interference and reflection effects was 
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Fig. 10. (a) Near-IR (left) and (b) mid-IR (right) transmission spectra and optical constants of crystalline C 2 H 2 . The vertical bar indicates the scale, and spectra are stacked for 
clarity. The samples giving the near- and mid-IR features shown had thicknesses of 4.6 and 0.49 pm, respectively. 


Table 2 

Positions, absorption coefficients (a), and band strengths (A) for selected C 2 H 2 features. 


Vibrational mode 3 

Amorphous C 2 H 2 at 12 K 



Crystalline C 2 H 2 at 1 5 K 



v (cm -1 ) 

2 (pm) 

oc (cm 1 ) 

A (10 18 cm molecule 1 ) 

v (cm -1 ) 

2 (pm) 

a (cm 1 ) 

A (10 -18 cm molecule -1 ) 

Vi + v 3 

6465 

1.523 

110 

0.42 

6442 

1.552 

87 

0.16 

5v 4 + 3v 5 

- 

- 

- 

- 

5174 

1.933 

82 

0.02 

Vi +v 5 

4072 

2.456 

129 

0.23 

4077 

2.453 

77 

0.11 

v 2 + 2v 4 + v 5 

3862 

2.589 

175 

0.21 

3864 

2.588 

115 

0.08 

v 2 + 2v 5 

- 

- 

- 

- 

3428 

2.917 

54 

0.11 

Vi 

3328 

3.005 

183 

0.59 

3333 

3.000 

57 

0.06 

v 3 

3240 

3.086 

22,470 

23.9 

3226 

3.100 

78,240 

35.0 

v 4 + v 5 

1371 

7.294 

1296 

3.1 

1393 

7.179 

5675 

6.7 

V 5 

742 

13.48 

11,740 

24.2 

770 

12.98 

19,460 

22.2 






761 

13.14 

14,780 



a See the footnote of Table 1 for comments on spectral assignments. 


agreement reached between the lab and the calculated results. 
Put another way, Fig. 11a corresponds to a calculation of an 
absorbance spectrum free of interference and reflection effects, 


while Fig. lib corresponds to an absorbance spectrum where such 
complications are present, just as they are in a laboratory ice 
sample. 


R.L. Hudson et all Icarus 228 (2014) 276-287 


285 


Table 3 

Band strengths for amorphous C 2 H 2 at 12 I< (10 -18 cm molecule -1 ). 


Vibrational mode 

A 

A' 

v 3 

23.9 

23.7 

v 4 + v 5 

3.1 

3.4 

V 5 

24.2 

27.5 


Table 4 

Band strengths for crystalline C 2 H 2 at ~70 K (10 -18 cm molecule -1 ). 

Vibrational A, Khanna 

mode et al. (1988) 

A, this 
work 

A', this 
work 

A', Khanna 
et al. (1988) 

A 1 , Dows 
(1966) 

v 3 

21.3 

35.8 

37.0 

35.7 

65.2 

v 4 + v 5 

5.98 

6.2 

7.6 

6.1 

8.47 

V 5 

13.8 

22.3 

26.9 

19.1 

18.3 




Wavenumber (cm' 1 ) 


Fig. 11 . Simulations of IR spectra of crystalline C 2 H 2 at 70 K as a function of 
thickness as calculated with (a) Beer’s Law, which uses only n(v) values, and (b) n(v) 
and /<(v) values to reproduce the interference and reflection effects in a real sample. 
Numbers in pm units refer to the ice’s thickness. 

Finally, we note that although the calculation of IR spectra 
from n(v) and k(v) values has long been studied, it is difficult 
to locate simple computer routines with sufficient detail and 
documentation for the necessary computations. We have found 
that the work of Swanepoel (1983) provides a particularly con- 
venient approach to calculating spectra from n(v) and k(v), and 
that while it has been heavily used by engineers and materials 
scientists, it seldom has been cited by workers within the astro- 
chemical community. A version of Swanepoel’s equations, 
adapted for Microsoft Excel (or equivalent) software, is posted 
on our web pages 1 for the use of all interested parties. Our n(v) 
and k(v) results, and those of others, can be used with this pro- 
gram to rapidly generate synthetic spectra without the need for 
more-expensive software. 


1 http://sdence.gsfc.nasa.gov/691/cosmicice/constants.html. 


5.3. Some comments on errors 

The input for each of our optical-constants calculations con- 
sisted of a high-quality IR transmission spectrum, a reference index 
of refraction (n 0 ), and the sample’s thickness (h). The good repro- 
ducibility of our results when data sets were combined from inde- 
pendent investigators working with different equipment over 
more than a year in two different laboratories gave us confidence 
that our IR spectra were of high quality and not a significant source 
of error. Our n 0 measurements resulted in small errors (~1%), and 
these can produce a vertical shift of the calculated infrared n values 
by the same amount, as suggested by Eq. (6). The thickness h of 
each ice sample could be measured with confidence to within 
about an eighth of an interference fringe, which corresponds to 
around 0.03 pm in thickness or about 6% of the thickness of the 
smallest samples (~0.5 pm). We believe that this is the dominant 
source of error in all of our results, as it is with those of Hudgins 
et al. (1993). To further check our error analysis, we also calculated 
optical constants starting with values of h and n 0 known to be 
incorrect by ~10%, and we then analyzed the resulting k values 
to determine band strengths (A) and absorption coefficients (a). 
In no cases were the resulting A and a of Tables 2-4 changed by 
as much as 10% even for acetylene’s strongest IR features in Table 2. 
Additional insight into uncertainties might be possible through a 
full propagation-of-errors analysis of the iterative Kramers-Kronig 
method that we used. So far as we are aware, no such analysis has 
appeared in the literature, despite the method’s long and wide- 
spread use. 

5.4. Applications 

The spectral properties of C 2 H 2 between 7000 and 500 cm -1 are 
of interest for comparison with both near-IR observations of cold 
icy surfaces, such as of TNOs, and mid-IR observations of cold icy 
interstellar grains. Our new data can be used to identify and quan- 
tify C 2 H 2 in ices, such as to estimate C 2 H 2 column densities from 
absorption spectra. Detailed discussions are available for the appli- 
cation of band strengths to astronomical spectra and the determi- 
nation of column densities (e.g., d’Hendecourt and Allamandola, 
1986). To take a specific case, spectral details and optical constants 
for acetylene’s v 5 band (742 cm -1 , 13.48 pm) are important since it 
is the strongest C 2 H 2 feature within the 5-20 pm observing range 
of the Spitzer-IRS instrument. Searches for IR spectral features of 
solid C 2 H 2 in dense cold molecular clouds are in progress, and by 
applying our laboratory data, information on the column density, 
phase, and possibly temperature history, of C 2 H 2 ices could be 
determined. 

The study of trans-neptunian objects also will be aided by our 
new data, both through remote sensing and spacecraft fly-bys. In 
particular, the Ralph/LEISA instrument on the New Horizons mis- 
sion will make hemispheric near-IR spectral observations at 
1.25-2.5 pm of Pluto and Charon at a resolution sufficient to 
map distributions of N 2 , CO, CH 4 , and possibly other trace spe- 
cies, such as C 2 H 6 and C 2 H 2 . As for TNOs other than Pluto, 
near-IR spectra of CH 4 and C 2 H 6 ices are consistent with the 
spectrum of Makemake, 2005 FY 9 (Brown et al., 2007) and possi- 
bly of Quaoar too (Schaller and Brown, 2007). Fig. 12 shows high- 
resolution H- and K - band spectra of Makemake (from Brown’s 
Figs. 5 and 6), where the thick line is a smoothed version of 
the data, and where the blue continuum and CH 4 ice have been 
subtracted (Brown, 2013, personal communication). 2 Fig. 12 also 
shows our laboratory spectra for amorphous C 2 H 2 , C 2 H 6 , and 


2 M. Brown has informed us that the modeling code used for the methane 
subtraction in our Fig. 12 differs slightly from that used in Brown et al. (2007). 
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C 3 H 8 (propane) at 30 K (Hudson et al., 2009 and unpublished 
work). As noted by Brown et al. (2007), their 100-pm C 2 H 6 grain 
model agrees well with Makemake’s spectrum beyond 2.25 pm 
and between 1.68 and 1.78 pm. Since the observed spectrum 
around 1.55 pm deviates from this same model (see Brown’s 
Fig. 6), extra absorbers are suggested, and C 2 H 2 is a candidate. In 
fact, acetylene is the only simple hydrocarbon that has a feature 
in this region. Makemake’s spectrum also deviates from the ethane 
model beyond 2.32 pm, and Brown et al. (2007) have suggested 
C 3 H 8 ice as a possible cause. Between 2.4 and 2.5 pm, C 2 H 2 , 
C 2 H 6 , and C 3 H 8 ices all have distinctive IR features, but telluric 
absorptions make TNO observations extremely noisy in this region 
(e.g., Brown’s Fig. 4). Of course for any quantitative spectral mod- 
eling, particle size and shape must be addressed and combined 
with spectra such as we show. 

Future observations from telescopes, such as SOFIA and NGST, 
above the Earth’s atmosphere may be needed to provide low-noise 
spectra. With such data we will have a better understanding of the 
contribution of C 2 H 2 and other hydrocarbons to the chemical 
inventory of icy surfaces of outer Solar System objects. 

6. Conclusions 

Here we have measured and presented, for the first time, spec- 
tra of amorphous and crystalline C 2 H 2 ice from 7000 to 500 cm -1 
(1.43-20 pm) at multiple temperatures in a single laboratory. 
From these data we have calculated the real (n) and imaginary 
(k) components of the complex index of refraction, from which 
measures of spectral intensity (A and a) can be computed. All re- 
sults are available on our web site 1 along with a spreadsheet pro- 
gram to calculate spectra. 
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